Dual quantum-correlation paradigms exhibit opposite statistical-mechanical properties 
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We report opposite statistical mechanical behaviors of the two major paradigms in which quantum 
correlation measures are defined, viz., the entanglement-separability paradigm and the information- 
theoretic one. We show this by considering the ergodic properties of such quantum correlation mea- 
sures in transverse quantum XY spin-| systems in low dimensions. While entanglement measures 
are ergodic in such models, the quantum correlation measures defined from an information-theoretic 
perspective can be nonergodic. 
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I. INTRODUCTION 

Quantum correlations - being one of the most impor- 
tant physical resources in quantum information science 
[1] - has led to the development of technologies based on 
quantum principles. There has been a continuous effort 
towards characterizing, quantifying, and utilizing quan- 
tum correlations for nearly the last two decades. Un- 
derstanding quantum correlations theoretically and their 
robust experimental generation remain important goals, 
although significant advances have been achieved in re- 
cent years on both fronts. 

Of the different measures that have been put for- 
ward to quantify and conceptualize quantum correla- 
tions, there are broadly two families. The first con- 
tains the ones that are defined within the entanglement- 
separability paradigm [2] , and they play a key role in de- 
scribing quantum information phenomena. Indeed, en- 
tanglement is the resource that drives the potentially 
revolutionizing applications in quantum information, in- 
cluding secure cryptography [3] , faster computers [4] , and 
better communication channels [5]. Recently, there has 
also been a considerable interest to study entanglement 
in quantum many-body systems [6-8]. 

Quantum correlations have also been conceptualized 
from an information-theoretical point of view - these 
forming the second family. Such correlations have the 
potential to be useful in understanding phenomena in 
which entanglement may not be necessary [9-11]. Impor- 
tant examples in this family are quantum discord (QD) 
[12] and quantum work-deficit [13]. Quite naturally, QD 
and quantum work-deficit have also been used to study 
many-body phenomena [14]. 

Investigations of ergodicity of a physical quantity in 
a many-body system plays an important role in under- 
standing statistical-mechanical properties of that sys- 
tem. The validity of a statistical-mechanical description 
of a quantity that characterizes a physical system de- 
pends on the behavior of that quantity in time-evolved 
states in comparison to that in the corresponding equi- 
librium states. A physical quantity is said to be er- 
godic if the time average of the quantity matches with 
its ensemble average. In this paper, we consider the 
status of ergodicity of quantum correlation measures in 
the quantum anisotropic transverse XY models, in low- 
dimensional spin-i lattices. The question of ergodicity of 



magnetization and classical correlations were considered 
in the quantum XY chain in Refs. [1-5-20]. Advances 
in the science of quantum information leads us to con- 
sider the question of ergodicity also for quantum correla- 
tion measures. Here we show that quantum correlation 
concepts defined from independent perspectives can ex- 
hibit opposing statistical-mechanical behaviors. We com- 
pare the ergodicity behaviors of quantum correlations 
defined within the entanglement-separability paradigm, 
such as concurrence [21] and logarithmic negativity [22], 
with those in the information-theoretic one, in the quan- 
tum XY spin systems in a transverse field in one di- 
mension, ladder, and two dimensions. We find that 
while entanglement measures display ergodicity irrespec- 
tive of the initial field strengths, quantum discord and 
quantum work-deficit are ergodic for both low and high 
transverse magnetic fields and nonergodic for moderate 
ones. For low fields, we prove ergodicity quite gener- 
ally: Arbitrary quantum correlation measures, irrespec- 
tive of whether defined in the entanglement-separability 
paradigm or by using information-theoretic concepts, are 
ergodic, for any interacting system in a non-commuting 
step function field, in any dimension. Moreover, we ana- 
lytically show that for high transverse magnetic fields, all 
quantum correlation measures are ergodic in the infinite 
quantum transverse XY chain for arbitrary anisotropy 
and arbitrary initial temperature. Exact diagonalization 
techniques are used to deal with high fields in higher 
dimensional systems. For moderate field strength, the 
Jordan- Wigner transformation is utilized in the case of 
an infinite one-dimensional (ID) system, while exact di- 
agonalization is used in higher dimensions. 

The paper is organized as follows. In Sec. II, we give 
a brief description of the model that we consider, and 
give simplified forms for the single- and two-site reduced 
density matrices of the model. In Sec. Ill, we define the 
logarithmic negativity and concurrence as measures of 
quantum entanglement and also define the information- 
theoretic measures, viz., quantum discord and quantum 
work-deficit. The notion of ergodicity is formally pre- 
sented in Sec. IV, and the same section also contains 
the general analytical ergodicity result for low fields. We 
study the ergodicity of the different quantum correlation 
measures for the infinite quantum spin chain in Sec. V, 
where we also present the analytical proposition for er- 
godicity at high fields. The ergodicity question for ID 
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finite spin cliains, of moderate sizes, is considered in Sec. 
VI, and we find tliat tlie answer is qualitatively the same 
as in the infinite case. Sections VII and VIII address the 
question of ergodicity, for the ladder and 2D quantum 
XY spin models respectively. Finally, we summarize our 
results in Sec. IX. 



II. THE XY MODEL IN A TRANSVERSE FIELD 

A. Description of the model 

An interacting Hamiltonian of spin-i particles on a 
lattice, with an external transverse magnetic field, is 
given hy H = Hiat — ^(O-^mag, where h{t) is a time- 
dependent transverse magnetic field. To ensure that the 
field part of Hamiltonian has a nontrivial effect on the 
system dynamics, one should have [i?int , -ffmag] 7^ 0. A 
simple way to obtain this is by choosing the field part of 
the Hamiltonian as iJmag ~ Tli^i^ ^"^^ interaction 
part as iJ^t = [(1 + 7)5f 5J + (1 - 7)Sf S-J] , 

with 7 7^ 0, where J is the coupling constant, 7 is the 
anisotropy measure, Sf , Sf , and are one-half of the 
Pauli spin matrices at the i*'^ lattice site, and (ij) in- 
dicates that the corresponding sum runs over nearest- 
neighbor pairs of the lattice. Periodic boundary condi- 
tions are considered in all the models dealt with in this 
paper. The total Hamiltonian studied, therefore, takes 
the form 

Hit) = J 5^ [(1 + i)s!sj + (1 - ^)sfsy] - hit) ^ SI 

(ij) « 

(1) 

This model is known as the anisotropic (quantum) XY 
model in a transverse field. The time-dependent mag- 
netic field h{t) is chosen as 

with a 7^ 0. Set hit) = h{t)/J, a~a/J. Although we 
carry out our investigations by considering a vanishing 
magnetic field for t > 0, the results remain qualitatively 
similar even for non-vanishing applied field strengths at 
nonzero time, i.e. for h{t) ^ for t > 0. 

B. Single- and two-site reduced density matrices 

In order to evaluate the quantities that are of interest 
in this paper, we need to find the single-site and two-site 
reduced density matrices of the equilibrium and evolved 
states of the system. 

A general single-site density matrix is given by 

where / denotes the identity on the qubit Hilbert space 
and M = tr[p^(T] is the magnetization vector, where a = 



(cr^, (T^, (T^) are the Pauli matrices. Let us first consider 
the single-site density matrix for the equilibrium state, 
Peq{t), at temperature T, where f3 = with ks being 
the Boltzmann constant. Due to symmetry, all the single- 
site density matrices of the equilibrium states are equal. 
We will denote them by p^git) (hiding the prefix /3). Now 
pi* it) = piqit), where the complex conjugation is taken 
in the computational basis, which (for each site) is the 
eigenbasis of the Pauli matrix . Therefore 

My^it) = tr[aVi,(t)] ^ 0. 

Moreover the Hamiltonian Hit) has the global phase flip 
symmetry i[H, Tli<7i] = 0), from which it follows that 

M^qit)^tr[a^piqit)]^Q. 

Consequently, the single-site density matrix of the equi- 
librium state reduces to 

pUt) = l{l + M^,it)a^). 

The evolved state does not necessarily have the property 
of being equal to its complex conjugation, and consider- 
ation of the global phase flip symmetry is complicated 
by the fact that the Hamiltonian is explicitly depen- 
dent on time. However, using the Wick's theorem, as 
in [17, 18, "Jo], the single-site density of the time-evolved 
state again turns out to be of the form 

p'it) = liI + APit)an. 

So, the single-site (transverse) magnetization of the equi- 
librium state is 

M4(0-tr[aV4(i)], 
while that for the evolved state is 

M'-it)=tT[a'p'it)]. 

The general two-site density matrix for the equilibrium 
as well as the time-evolved state can therefore be written 
as 

l[/®/+M^((7^®/+/®(T^)-f T'^a'^a')], 

where T*-' = tr[p^^(cr* ® a^)] represents two-site corre- 
lation functions. Once again by using Wick's theorem, 
we can show that the yz and xz correlations will van- 
ish for the evolved state, and for large time, also the xy 
correlation vanishes. For the equilibrium state, only the 
diagonal correlations, T^^, T^y , and T^^, remain. 

III. QUANTUM CORRELATION MEASURES 

A. Measures in entanglement-separability 
paradigm 

In this subsection, we will define two quantum cor- 
relation measures within the entanglement-separability 
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paradigm. Both have the virtue of being efScicntly com- 
putable for the case of two-qubit states, the latter being 
the quantum states at the focus of our study. 

1. Concurrence 

The concept of concurrence came from the definition 
of entanglement of formation and is proposed to quan- 
tify entanglement in two-qubit (possibly mixed) states 
[2 1] . The entanglement of formation of a bipartite quan- 
tum state is the amount of singlets, ";7|(|01) — |10)), that 
are required to prepare the state by local quantum op- 
erations and classical communication (LOCC). Here |0) 
and |1) form an orthonormal qubit basis. For two-qubit 
states, pab, there exists a closed form of the entangle- 
ment of formation in terms of the concurrence, which is 
given by G{pab) ~ max{0, Ai — A2 — A3 — A4}, where the 
Ai's are the square roots of the eigenvalues of pab Pab 
in decreasing order and pab = (cy <8) <^y)p*AB{^y ® 
with the complex conjugation being taken in the compu- 
tational basis. 



2. Logarithmic negativity 

Another measure of entanglement considered here is 
the logarithmic negativity (LN) [22]. The negativity, 
N[pab), of a bipartite state pab is defined as the abso- 
lute value of the sum of the negative eigenvalues of , 
where p^^ denotes the partial transpose of pab with re- 
spect to the A-part [24]. The logarithmic negativity is 
defined as 

EN{pAB)^\og^[2N{pAB) + l]- (3) 

For two-qubit states, p^"^ has at most one negative eigen- 
value [25]. Moreover, for two-qubit states, a positive LN 
implies that the state is entangled and distillable [24, 2G], 
while Em = implies that the state is separable [24]. 

B. Information-theoretic measures 

In this subsection, we will define two measures of quan- 
tum correlation defined from an information-theoretic 
perspective. In contrast to the entanglement measures 
defined in the preceding subsection, these are not com- 
putable in closed form. However, they can be efficiently 
computed, through numerical simulations, for two-qubit 
systems. 



1. Quantum discord 

Quantum discord is defined as the difference between 
two quantum information-theoretic quantities, whose 



classical counterparts are equivalent expressions for the 
classical mutual information [12]: 

Q{pab)^Apab)- J{pab)- (4) 

The total correlation, I(pab), of a bipartite state pab is 
given by [27] (see also [28, 29]) 

l^iPAB) = S{pa) + S{pb) - S{pab), (5) 

where S{g) = — tr(gilog2 g) is the von Neumann entropy 
of the quantum state g, and pA and p^ are the reduced 
density matrices of pab- On the other hand, JT{pab) 
can be interpreted as the amount of classical correlation 
in PAB, and is defined as 

J{pab) ^ S{pa) ~ S{pa\b)- (6) 

Here 

'S'(paib) = min^KS'(p^li), (7) 

is the conditional entropy of pab , conditioned on a mea- 
surement performed by B with a rank-one projection- 
valued measurement {Bi}, producing the states pA\i = 
^trB[(I^ (g) Bi)p{lA <8) Bi)], with probability p.^ = 
tr ab[{^a ® Bi)p{lA ® Bi)]. I is the identity operator on 
the Hilbert space of A. 

2. Quantum work-deficit 

Another important information-theoretic measure of 
quantum correlation is quantum work-deficit [13], which 
is defined for a bipartite quantum state as the difference 
between the amount of pure states that can be extracted 
under global operations and pure product states that can 
be extracted under local operations, in closed systems for 
which addition of the corresponding pure states are not 
allowed. 

The allowed class of global operations, known as 
"closed global operations" (CGO) are any sequence of 
(Gl) unitary operations, and (G2) dephasing of a given 
state PAB by using a set of projectors {Pi}, i.e., p — > 
Y.^P^PABP^, whcrc P^Pj = (5y P„ = I, with I be- 

ing the identity operator on the Hilbert space % on which 
Pab is defined. The number of pure qubits that can be 
extracted from pab by CGO is 

Ig{pab) = N - S{pab), 

where TV = log2(dim'H). The number of qubits that 
can be extracted from a bipartite quantum state pab 
under "closed local operations and classical communica- 
tion" (CLOCC), which consists of local unitary, local de- 
phasing, and sending dephased state from one party to 
another, is defined as 

Il{pab)^N- inf [S{p'a) + S{p'b)], (8) 
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where S{p'a) = tTB{A{pAB)) and S{p'b) = 

tYAiHpAB))- 

Quantum work-deficit is defined as the difference be- 
tween work, Ig{pab)i extracted by CGO, and that by 
CLOCC, Il{pab)- 

Mpab) = Ig{pab) - Il{pab)- (9) 

The quantity is not efficiently computable for arbitrary 
states, and therefore we restrict our attention to CLOCC 
consisting of measurement at a single party only. This 
restricted work-deficit is equivalent to quantum discord 
for bipartite states having maximally mixed marginals. 

IV. ERGODICITY 

In this section, we will formally describe the notion of 
ergodicity of a physical parameter, and will subsequently 
show that any quantum correlation measure is ergodic 
when the initial magnetic field is very low, irrespective of 
the other details of the Hamiltonian, and its dimension. 

To study the ergodicity of a physical quantity V in 
a given physical system, we evolve the system from an 
initial canonical equilibrium state at a given tempera- 
ture T, with the physical system being described by a 
Hamiltonian H = Hint — h(t)Hraag, where h{t) involves 
an initial kick to the system. The behavior of the time- 
averaged value, at large times, of "P, in the time-evolved 
state of the system is denoted by V°°{T, h{t), 7), where 7 
consists of the aggregate of all relevant parameters in the 
system, other than T and h(t). This time-averaged quan- 
tity is then compared to that of the canonical equilibrium 
state at large times and at temperature T', denoted by 
V'^'^T ,h{t = 00), 7). The physical quantity V will be 
called nonergodic, if for a given T, 

r°°{T,h{t),^) ^r'^''{T',h{t = 00),^) VT', (10) 

where T' is chosen in a physically relevant range. Other- 
wise, it will be termed ergodic. We term a physical quan- 
tity as strongly nonergodic if the relation (10) holds for all 
temperatures T'. Note here that the question of ergodic- 
ity in Ref. [19] was addressed by considering the effective 
temperature for a dynamical (time-evolved) state of the 
XY spin chain, which was determined by the condition 
that the equilibrium state must be on the same energy 
surface as that of the evolved state. Here we remove the 
energy constraint in investigating the ergodicity proper- 
ties of different physical quantities. In other words, the 
temperature T' is now not fixed by the energy (or any 
other) constraint, but is allowed to be within a certain 
relevant range of the initial temperature T. The length 
of this range can depend on several uncontrolled param- 
eters that may be active in an experimental realization 
of the theoretical model, including decoherence effects, 
disorder-induced effects, etc. 

The above definition of ergodicity is "microcanonical- 
like" , in the sense that for a particular physical quantity 



to be ergodic, one must have equality in Eq. (10). Such 
a definition is however of limited realistic applicability, 
as fluctuations of the corresponding physical values (on 
both sides of Eq. (10)) are disregarded. It is physically 
more meaningful to call a physical quantity as ergodic if 

(^-"(T', h{t = cx)), 7) - e, V--'\T', hit = 00), 7) + e), 

(11) 

for some T', and where the choice of e is to be elaborated 
below. A physical quantity which is ergodic according 
to Eq. (11) will be termed as "canonically ergodic". A 
quantity that is ergodic (i.e. it satisfies Eq. (10) with 
an equality, for some T') is of course canonically so. The 
choice of e will depend on the physical system under con- 
sideration, and the physical quantity that we are study- 
ing. Typically, it will be of the same magnitude as the 
error in evaluating that quantity in an experimental im- 
plementation of the system. Where such an error is not 
apparent, one may choose e as the standard deviation 
of the physical quantity in the ensemble that is being 
considered. 

Proposition 1. Consider the Hamiltonian H = Hint — 
h{t)H„iag, where [Hint.Hmag] 7^ 0. The time-dependent 
"magnetic field", h{t), is a step function and is of the 
form h{t) = a,t <0 and h{t) = 0,t > 0. For low initial 
field a, i.e. in the neighborhood of a « 0, all quantum 
correlation measures are ergodic, in any dimension. 
Proof. If we choose a to be zero which is the same as the 
magnetic field for t > 0, then there will be no evolution 
in the system and hence all the quantum correlations of 
the evolved state will be the same as that of the initial 
thermal equilibrium state. In the vicinity of a « 0, and 
with the initial temperature T, it is always possible to 
find a temperature T « T for which all quantum corre- 
lation measures of the evolved state will be same as that 
of the thermal equilibrium state at T . Here we assume 
that Hint and H^ag do not depend on time. □ 
Although we state the proposition for low fields, the 
theorem holds for any value of a which is close to a con- 
stant h{t) for i > 0. 



V. INFINITE QUANTUM XY SPIN CHAIN IN 
A TRANSVERSE FIELD 

We now study and compare the statistical-mechanical 
behavior of quantum correlation concepts defined in 
the two qualitatively different paradigms, viz., the 
entanglement-separability and the information-theoretic 
ones, as described in Sec. HI, for the infinite ID XY 
model (see Eq. (1)). From Prop. 1, we obtain that for 
low fields, the measures from both the paradigms will 
exhibit ergodicity. Here we show, analytically, that for 
high fields, they are also ergodic. In contrast, moder- 
ate magnetic fields produce opposite effects in the two 
paradigms, as we will show later. 
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Proposition 2. For high fields i.e. for a oo, any 

quantum correlation measure is ergodic for the ID trans- 
verse XY model for arbitrary initial temperature and the 
anisotropy. 

Proof. The ID transverse XY model lends itself to ex- 
act diagonalization [17, 18, 23]. We perform the cor- 
responding Jordan- Wigner and Fourier transformations 
successively and find that the transverse magnetization 
in the time-evolved state at long times, considering that 
the initial state is the canonical equilibrium state at tem- 
perature r, is given by 



1 r ta.rh(W)/2) 
ttJo ^ A(5)A2(0) 
X cos 0[(cos 4> ~ a) cos + 7^ sin^ 



(12) 



Here K{x) = {7-^ sin^ + [a; — c 
The correlations are given as follows: 



12}% and P = /3J. 



rj^xx ^ G{-i), ryy = g(i), t'' = []VPf~G{i)G{~i) 

where G{R) (for R = ±1) is given by 

I r tanh(/3A(a)/2) , , , , , 

TT Jo A(a)A2(0) 

x(7^ sin^ (/) + (cos(/) - a) cos(^). (13) 

Moreover, T'^y = T^'^ = 0. As a 00, G(i?) for 
i? = ±1 for all 7 and ^. Therefore, also T^^ ~ Tyy = 0, 
while T^z reduces to 4[M'']2. And 



I + I7I 



as a — > 00. Therefore for high fields, the two-site reduced 
density matrix turns out to be of the form pA® Pb, with 
PA = PB = j{I + M^a^), and hence all quantum corre- 
lation measures of the evolved state vanish. Vanishing of 
quantum correlation measures is also true for the ther- 
mal equilibrium state for T — >■ cxd, implying an ergodic 
nature for all the measures of quantum correlations. It 
also shows that is nonergodic for any choice of 7, if 
one chooses h{t) = for i > 0. □ 



Ergodicity of entanglement measures 

Combining Propositions 1 and 2, we find that all quan- 
tum correlation measures behave similarly for low and 
high fields. The situation is drastically different for mod- 
erate fields. Wc begin by demonstrating that two en- 
tanglement measures, viz., concurrence and logarithmic 
negativity, exhibit ergodic behavior, also for moderate 
fields. Let us begin by considering the behavior for con- 
currence (see Fig. 1 (top-left)). The computations are 
performed by using Jordan- Wigner transformations. For 
the purpose of presentation in all the figures in this paper, 
we choose the anisotropy 7 to have the value ^, and the 




FIG. 1. (Color online.) Ergodic entanglement measures ver- 
sus nonergodic information-theoretic measures. Entangle- 
ment measures (concurrence and logarithmic negativity, in 
ebits) and information-theoretic quantum correlation mea- 
sures (quantum discord and quantum work-deficit, in bits) 
are plotted on the vertical axes against the dimensionless pa- 
rameter j3, for the infinite quantum transverse XY spin chain. 
Concurrence (C, top-left) and logarithmic negativity (-Ejv, 
top-right) are ergodic for any value of initial transverse field 
a. In contrast, information-theoretic quantum correlations - 
quantum discord [Q, bottom-left) and quantum work-deficit 
(A, bottom-left) - are nonergodic for moderate values of the 
applied transverse field. The curves with black circles are for 
the canonical equilibrium state at large time, as functions of 
j3. We have plotted the values of the corresponding physical 
quantity, for the time-evolved state with the initial state of 
the evolution being the canonical state at /3 = 20, as hori- 
zontal lines for different values of the applied transverse field. 
Precisely, the applied fields used for the horizontal lines are 
a/J = 0.2 (red squares), a/J = 0.6 (blue diamonds), and 
a/J = 2 (pink triangles). Note that a/J is a dimensionless 
quantity. Here 7 = 0.5. 



temperature of the initial (canonical) state of the time- 
evolution to be such that /? = 20. These choices are of no 
particular significance, and all the results obtained have 
been confirmed to be qualitatively similar for all nonzero 
values of 7 and for all low initial temperatures. A nonzero 
anisotropy is necessary for the change in the transverse 
field to have any nontrivial effect on the evolution, and 
high temperature canonical states are known to have zero 
or near-zero quantum correlations. For the evolved state, 
the time-averaged values of concurrence, at large times, 
for the exemplary values of a/J=0.2, 0.6, and 2.0 are 
0.153, 0.1149, and respectively. These are depicted in 
the figure (Fig. 1 (top-left)) as horizontal lines at the 
corresponding heights. The concurrence of the canoni- 
cal equilibrium state starts off from zero for /? = and 
converges to 0.153 at low temperatures. Hence, there 
always exists a temperature for which the correspond- 
ing equilibrium value matches the time-averaged value 
in the evolved state, up to the 3rd significant digits. 
The feature is shared by other values of a/J as well. 
Therefore we conclude that the concurrence is ergodic 
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for the infinite ID XY spin model. For low values of 
the initial field, the concurrence may only be canoni- 
cally ergodic. It is clear from Fig. 1 (top-right) that, 
along with the concurrence, the logarithmic negativity, 
another quantum correlation measure defined within the 
entanglement-separability paradigm, is also ergodic (or 
at least canonically so) in this model. 

In Ref. [19], the question of ergodicity of logarithmic 
negativity was addressed for the ID XY spin model by 
requiring that the energy of the equilibrium state and the 
dynamically evolved state is the same. This, however, is 
a stringent condition, and is not usually considered to be 
necessary for ergodicity considerations (see e.g. [17, LS]). 
Rather surprisingly, we find that by removing this en- 
ergy constraint, entanglement-like quantities (including 
logarithmic negativity) become ergodic in this model. 

Nonergodicity of information-theoretic measures 

Measures of quantum correlation defined from an 
information-theoretic perspective exhibit a behavior that 
is opposite to that of entanglement measures, for mod- 
erate fields. Let us first consider the case of QD. For 
low values of the applied transverse field, the QD re- 
mains ergodic (see Fig. 1 (bottom-left)). This is just 
like for entanglement measures and also follow directly 
from Prop. 1. However as the applied transverse field 
is increased to moderate values, there is no temperature 
for which the time-averaged value of QD matches with 
that of the canonical equilibrium value - QD is strongly 
nonergodic in the vicinity of a = 0.6. Again the compu- 
tations are performed by using successive Jordan- Wigner 
and Fourier transformations. For even higher magnetic 
fields - e.g. for a/ J = 2.0, - QD attains ergodicity only if 
we allow temperatures much higher - precisely two orders 
of magnitude higher - than the initial temperature, and 
hence we still term the situation as nonergodic, although 
not strongly so. 

Surprisingly, we have also observed that QD shows 
nonergodicity even when its constituent quantities, I and 
JT", are both ergodic. This arises due to the fact that 1 
and S become ergodic at different temperature ranges, 
and so the ergodicity of QD does not remain a linear 
function of the ergodicities of 2 and J. QD must be 
ergodic for a — > oo, follows from Prop. 2. Very similar 
features are observed in the other information-theoretic 
quantum correlation measure, viz., quantum work-deficit 
(see Fig. 1 (bottom-right)). 

Discussion 

In an attempt towards understanding the opposite be- 
havior of the quantum correlations in the two prominent 
paradigms, let us begin by contrasting QD against loga- 
rithmic negativity for the evolved state, with respect to 
the initial transverse field at large time (see Fig. 2). We 
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FIG. 2. (Color online.) Quantum discord versus logarithmic 
negativity at large time. Quantum discord (in bits, left figure) 
and logarithmic negativity (in ebits, right figure) of the time- 
evolved state at t — >■ oo are plotted on the vertical axis, for the 
infinite quantum transverse XY spin chain is plotted against 
the dimensionless transverse magnetic field parameter a/ J , on 
the horizontal axis. The temperature of the initial canonical 
state of the evolution is taken such that /3 = 20, and 7 = 0.5. 



find that while logarithmic negativity is monotonic (de- 
creasing) for the whole range of a, QD is increasing for 
low values of magnetic field and decreasing for high val- 
ues. This different behavior of the two measures hints at 
the different statistical-mechanical behavior for QD and 
logarithmic negativity. 

Another way to interpret the opposite behavior ob- 
served is the following. Classical correlators, given by 
Tij = {ai ® Uj), of the states of a many-body quantum 
system has long been known in the literature of con- 
densed matter [.:>()] and Bell inequalities [31], among oth- 
ers, as important characterizing quantities of the phys- 
ical system. Quantum correlations have recently been 
identified to be interesting quantities of many-body sys- 
tems [0-S]. A typical multisite quantum state contains 
both quantum and classical correlations. The measures 
of entanglement strives to filter out only the quantum 
part of the correlations in a way that is very stringent, 
as compared to the information-theoretic measures of 
quantum correlations. Indeed, even separable states have 
nonzero quantum discord [12] and quantum work-deficit 
[13]. Therefore, it is reasonable that the information- 
theoretic measures of quantum correlations will behave 
more closely like the classical correlators, while measures 
defined within the entanglement-separability paradigm 
will behave oppositely thereof. The classical correlators 
in these systems have long been known to be generally 
nonergodic [15-18] (cf. [20]). They are ergodic for low 
transverse fields, while becoming nonergodic for mod- 
erate fields - exactly the feature obtained for quantum 
discord and quantum work-deficit, while exactly the op- 
posite to that obtained for concurrence and logarithmic 
negativity. 

Magnetic field regions.- It is clear from Fig. 2 that 
whenever the applied field strength a/ J is approximately 
less than 0.4, all the quantum correlation measures, 
whether they belong to the entanglement-separability 
paradigm or the information-theoretic one, have a low 
derivative (slow change) with respect to the applied field. 
We call such regions as "low field regions". Similar be- 
havior (of the measures) is also observed when a/ J is ap- 
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proximately above 1.0. Such regions are called the "high 
field regions" . The remaining fields (on the positive field 
axis) fall in what we term as the "moderate field region" . 
In this moderate field region, all the quantum correlation 
measures have a substantially larger derivative values (in 
comparison to the same in the other regions) with respect 
to the applied field. It may be noted here that the quan- 
tity a/ J is a dimcnsionless quantity, so that the field (a) 
regions arc already defined in energy (J) units. 



VI. FINITE QUANTUM XY SPIN CHAIN IN A 
TRANSVERSE FIELD 

In the previous section, we considered the infinite 
XY spin chain in a transverse field and found that 
information-theoretic quantum correlation measures and 
those of the entanglement-separability paradigm have 
countering statistical-mechanical behavior. It will be in- 
teresting to see whether such complementary behavior is 
sustained even in higher-dimensional systems. However, 
the analytical methods do not exist for a ladder or a two- 
dimensional model and so we have to use numerical sim- 
ulations to study them. The question remains whether 
results obtained by exact diagonalization can infer cor- 
rectly the ergodic behavior of ladder or two-dimensions. 
Hence, we first consider the finite XY spin chain of mod- 
erate length and compare their results with the infinite 
one. 



C El 




FIG. 3. (Color online.) Reproducing the content of Fig. 1 
for finite chain of 12 spins with periodic boundary conditions. 
All considerations remain the same as in Fig. 1. 

The Hamiltonian for a spin chain of 12 spins, with 
periodic boundary conditions, is given by 

12 12 

h = jY,[{i + i)s^st+^ + (1 - i)sfsy^,] -h{t) SI 

i=l 1=1 

(14) 

Similar to the case of the infinite spin chain, an initial 
transverse field is applied as a disturbance on the initial 



canonical equilibrium state at /? = 20, and then we al- 
low the system to evolve. Later, we compare the long 
time-average of the physical quantities for the evolved 
state and the canonical equilibrium state at different 
field strengths. The plots of concurrence, logarithmic 
negativity, quantum discord, and quantum work-deficit 
for the evolved state and the canonical equilibrium as 
functions of temperature are given in Fig. 3. Simi- 
lar to the case of infinite XY chain, we observe that 
entanglement-based measures, viz., concurrence and log- 
arithmic negativity, remain ergodic for the entire range 
of the applied initial transverse magnetic field. However, 
the information-theoretic quantum correlation measures 
are ergodic for low and high values of applied transverse 
field and nonergodic at moderate values. To illustrate 
this, let us consider the case of QD, which remain ergodic 
for low initial fields (Fig. 3 (bottom-left)). For moder- 
ate values of the initial field, QD becomes nonergodic, 
i.e, there is no temperature for which the time-averaged 
value of QD for the evolved state is equal to the QD of 
the equilibrium state. Hence, the statistical-mechanical 
behavior of the quantum correlation measures defined 
in the entanglement-separability paradigm (concurrence 
and logarithmic negativity) and that of the measures de- 
fined in the information-theoretic one (quantum discord 
and quantum work-deficit) for finite ID XY model of 
length N = 12, mimics the properties of the infinite 
spin chain. This induces us to study such statistical- 
mechanical behavior in the case of similar Hamiltonians 
in higher dimensions, where we will be forced to consider 
finite-sized systems, as these systems are not analytically 
diagonalizable. This will be carried out in the two suc- 
ceeding sections. 



VII. QUANTUM XY LADDER MODEL IN A 
TRANSVERSE FIELD 

The ladder arrangement of spins can be thought of as 
an intermediate dimension between one and two dimen- 
sions. These systems are interesting as they have revealed 
several intriguing features of both one and two dimen- 
sions. Ladder models have been successfully realized in 
laboratories in certain materials and currently available 
techniques in cold gas systems also offer the possibility of 
realizing such models [32]. The ladder model with peri- 
odic boundary conditions and XY interactions on steps 
as well as rungs is considered. 

Again the system is evolved from the canonical equi- 
librium state with an initial temperature oi P — 20, after 
applying a kick in the transverse magnetic field. As in 
the case of the infinite and finite spin chains, for mod- 
erately high fields, QD and quantum work-deficit are 
nonergodic, while nearest-neighbor entanglement (con- 
currence as well as logarithmic negativity) of steps as 
well as rungs remains ergodic. From Fig. 4 (right), it is 
clear that the evolved state logarithmic negativity, plot- 
ted as horizontal lines corresponding to transverse field 
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FIG. 4. (Color online.) Ergodicity versus nonergodicity of 
quantum correlation measures on XY ladders. Quantum dis- 
cord (left, in bits) and logarithmic negativity (right, in ebits) 
for a ladder of eight spins with periodic boundary condition. 
The considerations are the same as in Fig. 1, except for the 
orange upside-down triangles which are for a/ J — 0.1. 



strengths of a/J=0.1, 0.6, and 2.0, wiW definitely cross 
the logarithmic negativity of the canonical equilibrium 
state for some values of the temperature. Such crossing 
of the equilibrium and evolved state entanglement curves 
does not happen for QD for initial fields a/ J > 0.15 
(Fig. 4 (left)). Thus the statistical behavior of nearest- 
neighbor entanglement-based measures and information- 
theoretic ones, qualitatively matches with those of the 
XY spin chain: entanglement measures are ergodic while 
information-theoretic ones are not so, for moderate field 
strengths. 



VIII. TRANSVERSE 2D XY SPIN CHAIN 

The two-dimensional quantum XY model with a trans- 
verse magnetic field plays an important role in both 
many-body physics and quantum information. In par- 
ticular, the dynamics of this model in the case of unit 
anisotropy is used for realizing the measurement-based 
quantum computation [3.3]. Here we consider a two- 
dimensional array of 12 spin-i particles arranged on a 
square lattice, and interacting via nearest-neighbor XY 
interactions with periodic boundary conditions. The ar- 
rangement therefore forms a torus. 




measures show ergodicity also in this case. However, we 
find that for a/ J = 0.6, the time-averaged QD of the 
evolved state has the value 0.0433 for an initial tem- 
perature, corresponding to P = 20, while the QD of 
the canonical equilibrium state converges to 0.0424, and 
hence leads to nonergodicity. The same feature is shared 
by quantum work-deficit. In contrast, entanglement- 
based measures remain ergodic for all values of the trans- 
verse field, including a/ J = 0.6. 

One may expect that the ergodic behavior of the quan- 
tum correlations in the 2D model will be different from 
that in the ID model. However, simulations in finite 2D 
systems show a qualitatively similar behaviour as in ID. 
However, some classical correlations (of the form T*-') 
turn out to be ergodic in 2D, while they are not so in ID 
[20]. 



IX. CONCLUSION 

We have considered the statistical-mechanical proper- 
ties of quantum correlations of the transverse quantum 
XY model with nearest-neighbor interactions in low di- 
mensions with periodic boundary conditions. We have 
shown that any quantum correlation measure of a Hamil- 
tonian with a time-independent interaction and a non- 
commuting time dependent field is ergodic for low fields, 
irrespective of the other parameters of the Hamiltonian. 
Moreover, for the infinite XY spin chain, with a high 
initial transverse magnetic field, we have analytically 
demonstrated that any quantum correlation measure is 
again ergodic. Furthermore, we found that the quan- 
tum correlation measures defined in the entanglement- 
separability paradigm, equilibrate with time, and hence 
remain ergodic, while the information-theoretic quan- 
tum correlation measures can be nonergodic for mod- 
erate initial fields. The thesis is true for the infinite 
chain, which is exactly diagonalizable. Low-dimensional 
finite systems, that cannot be analytically diagonalized - 
ladder and two-dimensional models with periodic bound- 
ary conditions - mimics the one-dimensional infinite one. 
We therefore obtain a distinct partition that divides 
the world of quantum correlation measures - ones that 
are ergodic as opposed to the ones that are not so, 
in low-dimensional quantum transverse XY spin-i sys- 
tems. More interestingly, the partition coincides with 
that which divides them into entanglement-based quan- 
tum correlation measures and information-theoretic ones. 



FIG. 5. (Color online.) Ergodicity versus nonergodicity on 2D 
XY models. Quantum discord (left, in bits) and logarithmic 
negativity (right, in ebits) for a 2D array of 12 spins with 
periodic boundary condition. All considerations are the same 
as in Fig. 1. 

For various transverse field strengths, viz., a/ J ~ 
0.2,0.6, and 2.0, we find that the entanglement-based 
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